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ABSTRACT 

Analysis of a database of solar coronal mass ejections (CMEs) and associated flares over the period 
1996-2007 finds well-behaved power law relationships between the 1-8 A flare X-ray fluence and CME 
mass and kinetic energy. We extrapolate these relationships to lower and higher flare energies to 
estimate the mass and energy loss due to CMEs from stellar coronae, assuming that the observed X-ray 
emission of the latter is dominated by flares with a frequency as a function of energy dn/dE = kE~°' . 
For solar-like stars at saturated levels of X-ray activity, the implied losses depend fairly weakly on the 
assumed value of a and are very large: M ~ 5 x 1O~^°M0 yr~^ and E ~ Q.ILq. In order to avoid such 
large energy requirements, either the relationships between CME mass and speed and flare energy 
must flatten for X-ray fluence > 10'^^ erg, or the flare-CME association must drop significantly below 
1 for more energetic events. If active coronae are dominated by flares, then the total coronal energy 
budget is likely to be up to an order of magnitude larger than the canonical lO^^Lf,,,; X-ray saturation 
threshold. This raises the question of what is the maximum energy a magnetic dynamo can extract 
from a star? For an energy budget of 1% of Lboi, the CME mass loss rate is about 5 x IQ^^'^Mq yr^^. 
Subject headings: Stars: winds, outflows — Sun: coronal mass ejections (CMEs) — stars: flares — 
X-rays: stars 



1. INTRODUCTION 

The rate of mass loss from unevolved late-type stars 
is notoriously difficult to constrain. In case of the 
Sun, the wind can be directly observed and sampled by 
spacecraft, and amounts to a mass loss rate of about 
2 X lO^^^Af0 yr"-'^. Such a weak flow of ionized gas 
from other stars cannot be detected directly using in- 
strumentation available today, and there are presently no 
direct detections or measurements of winds from solar- 
like stars. Mass loss rate upper limits based on radio 
obse rvations are in the ra nge of sever a l \G~ ^^Mp, yr~^ 
(e.g. iGaidos et al1[2000l ). iWood et al.l (|200l devised a 
method of indirect assessment based on H Lya absorp- 
tion due to interstellar H I that is heated in the inter- 
action region between the wind and the local interstel- 
lar medium. They estimated rates in the range 10~^^- 
1O~^^M0 yr~^, with evidence for higher mass loss rates 
for more X-ray luminous stars. 

Theoretical progress in predicting solar-like winds has 
been hampered by a persistent lack of understand- 
ing of the basic mechanisms responsible for producing 
them. Inspired by the recent success of turbulence- 
driv en coronal heating and s olar wind acceleration the- 
ory, iCranmer fc Saaii ()201lD developed a wind model 
based on the energy flux of magnetohydrodynamic tur- 
bulence from the subsurface convection zone. For a solar- 
like star, they predict mass loss rates that decline steadily 
from a few IQ^^'^Mq yr^^ at the zero-age main-sequence 
and rotation periods of 1 day to ~ 1O~^^M0 yr~^ at 
rotation periods of 60 days or so, in re a sonab le agree- 
ment with the estimates of iWood et al.l ()2002l see also 
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IWood et al.l[2005l) . 

Further progress in understanding mass loss of main- 
sequence late-type stars is strongly motivated by 
the effect winds have on stellar rotation evolution 
(e.g. ^cbcr fc Davis 1967; Stauffcr fc Hartmann 19&_ 
Kawalei<1988tlMatt fc Pudritzi2008HReiners fc Mohant 



20llr and consequently on stellar magnetic activity (e.g. 
Pallavicini et al. | Il981t IWright et all | 201lD . and on in- 



terplaneta ry medium env ironments (e.g. iPreusse et al.l 
2005; Lam mer et 311120071 ) . 

One aspect of mass loss that remains to be thor- 
oughly investigated on other late-type stars is that 
due to coronal mass ejections (CMEs). On the Sun, 
CMEs are observed to eject from 10^^ to lO^^g of 
ma gnetized plasma into the interplanetary rnedium 
(e.g|YMhirQ_&Gopalswamv''2009t iVourhdas et~aIll2010D . 
The integrated mass loss from CMEs ca n amount to sev- 
eral p ercent of the steady wind rate (e.g lVourlidas et all 
|2010( ). At first sight this suggests CMEs are going to 
be of little importance in the stellar context. However, 
on the Sun CMEs are associated with fiares, and mag- 
netically active stars are widely int e rpreted to b e dom- 
inated by flares (e.g. iGuedell Il997l: iDrake et al.ll2000D . 
Since the most magnetically active solar-like stars can 
attain X-ray luminosities more than 1000 times that of 
the Sun, there is scope for vigorous CME activity, espe- 
cially in the context of recent giant flare detections on 
solar-type stars based on optical Kepler photometry by 
iMaehara erahl (|2012D . The importance of CMEs on ac- 
tive stars has also been raised in the context of erosion of 
the atmospheres of "Hot Jupiters" (jKhodachenko et al.l 
jp07a) and the habi tability of planets around M dwarfs 
(fkhodachenko et al.l 120079 ). while lAarnio et al.l ()2012[ ) 
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suggest that CMEs associated with flares could be an 
important contribution to angular momentum loss on 
pre-main sequence stars. 

Here, we examine the consequences for mass loss and 
energy loss of flare-dominated coronae based on extrapo- 
lation of the observed behavior of a large sam ple of CMEs 
compiled bv lYashiro fc Gopalswanivl ()2009f ). The impli- 
cations are quite striking and provide an indirect means 
to begin to assess how CMEs might behave on stars much 
more magnetically-active than the Sun. 

2. MASS AND ENERGY OF SOLAR CORONAL 
MASS EJECTIONS 

As a guide to t he CME behaviour of ac t ive st ars we 
look to the Sun. lYashiro fc Gopal swamvl ([2009') stud- 
ied the statistical relationships between solar flares and 
CMEs observed over the period 1996-2007 (see also 
iVourlidas et al.l 120101 for a description of some of the 
observational aspects of CME measurements). They 
compiled a database of soft X-ray flares observed by 
the Geostationary Operational Environmental Satellite 
(GOES) that were associated with CMEs observed by the 
Large Angle and Spectrometric Coronagraph (LASCO) 
on board the Solar and Heliospheric Observatory (SoHO) 
mission. The CME-flare association fraction was ob- 
served to increase with flare peak X-ray flux, fluence, 
and duration, as had been noted in earlier studies (e.g. 
I Andrews! I2003f ). and a good correlation was found be- 
tween the flare fluence and the CME kinetic energy. 

The distribution of CME ejected masses as a function 
of t he associated flare GOES 1-8 A X-ray fluence from 
the lYashiro fc Gopalswamvl ()2009f ) sample is illustrated 
in Figure [1] While the data show a very large scatter 
in CME mass at a given flare energy, we find that the 
mean of these data over small X-ray fluence bins is well- 
behaved and adheres quite closely to a power law. Such 
a power law relation between ejected mass and the peak 
X-ray flux of the as s ociate d flare has also been pointed 
out bv lAarnio et al.l (|2011l ). A linear fit to the logarithm 
of the variance- weighted means of 20 point bins yields the 
following power law relationship (in cgs units) between 
ejected mass and fiare fluence, 

m,{E) = ^lE^ ■ 

10-l•5T°•^ /3 = 0.59 ±0.02. (1) 

The constant of proportionality and power law index 
uncertainties are strongly anti-correlated. They were 
determined using a M onte Carlo multiple imputation 
bootstrap ()Rubin|[l996[ ) in which the distribution of fit 
parameters was estimated from repeated re-fitting of a 
randomly-drawn 2/3 of the data augmented to the full 
s ample size by random dr a ws fro m this sub-sample. 

Yashiro fc Gopal swamvl ((2009) found a similar distri- 
bution of CME kinetic energy as a function of associated 
flare X-ray fluence. The data are shown in Figure [TJ to- 
gether with the histogram of the mean of every 20 points 
and a power law fit to this. From this fit we find the mean 
CME kinetic energy to vary with flare X-ray fluence as 
(again in cgs units) 

Eke =- vE^ ; 

= 100-8i=F0.85^ ^ = 1.05 ±0.03. (2) 
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Fig. 1. — CME mass (top) and kinetic energy (bottom) vs X- 
ray flu ence of the associated flare from the Yashiro_45_ Gopalswair^ 
l|2009l ) sample. The green histograms are the means over twenty 
data points and the red lines are linear flts to these means. The 
dashed blue lines are the linear flts multiplied by the CME-flare 
association rate given by Eqn. [3] In the upper panel, the dashed 
grey line follows a constant ratio of mass loss to GOES X-ray en- 
ergy loss, expressed as rates, M = 10~"'^''(Lx/10^'')Mq yr~^. In 
the lower panel, the dashed grey line represeniis equivalence of the 
kinetic and X-ray energies. The red line in this panel corresponds 
very closely to a factor of two hundred times the X-ray fluence. 



Again, the uncertainties were determined using a Monte 
Carlo multiple imputation bootstrap and are essentially 
anti-corrclatcd. Also shown in Figure [T] is the locus of 
equivalence between X-ray and kinetic energies. The lat- 
ter lies about a factor of 200 above the former, indicating 
that, for a given flare, the energy release will be totally 
dominated by the energy of the associated mass ejection. 
We return to the consequences of this in ^ 

We will find it useful below to also express the CME- 
fiare association fractio n as a power law. We find the 
lYashiro fc Gopalswamvl ()2009[ ) association fraction as a 
function of X-ray fluence, f{E), can be well-represented 
by 

f{E) = 1 for E >3.5x 10^^ erg 
f{E) = C^^'^ for E <3.5x lO^^erg; 

C = 7.9 X 10"^^ 6 = 0.37. (3) 
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3. ESTIMATING TOTAL STELLAR 
CME-ASSOCIATED ENERGY AND MASS LOSS 

Flare occurrence in the coronae of the Sun and stars 
has been shown by a number of studies to follow a power 
law distribution in frequency as a function of flare energy 
of the form 



where fc is a normalization constant (e.g. 'Drake' '19711; 
Datlowe et all 1974; Lin et al. 1984; Hudson 1991; 
Bai 1993; Porter et al. 1995; Krucker & Bcnz 199& 
Audard et al.„ ,2000: Ka shvap et al.. .2002: Giide l et all 
20031 : iHannah et al.i[20l'i'b . Analyses of stellar EUV and 



X-ray light curves furthermore suggest active stellar coro- 
nae are dominated by a superposition of flares, and tend 
to find values of the frequency vs. energy power law in- 
dex in the range a = 2-2.5 for all stellar types, includ- 
ing dwarfs with spec t ral type G-M ([A udard et al. 2000; 
Kashvap et all [200l iGiidel et al .' '2003; Telleschi et al.i 



20051 see also the earlier work of CoUu ra et al 
and T Tauri stars (jCaramazza et al.ii2007t Istelzer et al.l 
|2007[ ). Similar indices were o bserved for opti c al flar es de- 
tected on solar- type stars bv iMaehara et all (|2012[ ) from 
Kepler photometry. 

The total flare power is given by the integral over min- 
imum and maximum flare energies 



P = 



EkE-" dE = 



2- a 



[El 



2-a 

max 



Et 



(5) 



Taking the X-ray luminosity, Lx , as the observable proxy 
for the flare power, the constant k is given by 



Lx{2-a) 



(6) 



For the case in which CME mass loss is a function of 
the associated flare X-ray fluence, the total mass loss rate 
from CMEs is 



E„ 



m,[E)f{E)^dE. (7) 



When combined with Eqns.[Tl[3]and|ni the resulting rate 
is 



M, = ^JLC,Lx 



1 + (3 + 5 -a 



En 



Et 



which we can now evaluate for suitable choices of the 
minimum and maximum flare energies. Similarly, the 
total CME-associated kinetic energy loss rate is 

•m n T 'm in 



Eke = vCLx 



a 



(9) 

Mass loss and kinetic energy loss rates are illustrated 
as a function of the power law index a in Figure [5] for 
different values of Emin and Emax- Here, we have nor- 
malised to a coronal luminosity of Lx — 10'^° erg s~^ 
in the GOES 1-8 A bandpass, which, for a fairly typical 
coronal temperature for the most active solar-like stars 
of 2 X 10'' K, corresponds to Lx ~ 3 x 10^*^ erg s~^ in 
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Fig. 2. — CME mass (top) and kinetic energy (bottom) loss rates 
vs. power law index a for a 1-8 A X-ray luminosity of Lx = 
10^0 gj-g g-i^ according to Eqns. |8] and |9l respectively. Note that 
there are singularities corresponding to power law indices a = I + 
(5 + 5 and a = 2. The grey shaded areas represent the uncertainties 
in the loss rates corresponding to the uncertainties in the power law 
fits in Eqns. [T] and [2] 



the 0.2-2.5 keV and 0.5-10 keV bandpasses. This X- 
ray luminosity is that of a coronally-saturated solar-like 
star with a ratio of X-ray to bolometric luminosity of 
Lx/Lboi ~ 10~'^i suc h as 47 Gas B or EK Dra (see, e.g., 
iTelleschi et"aI|[20Q5l ). 

For fiducial limits we have adopted Emax — 10'^^ erg 
and Emin — W^^Emax- The former corresponds to a 
reasonably large but fairly common flare on an active 
solar-type star. Figure [2] demonstrates that the particu- 
lar choice of these integration limits is not important — 
changing the lower limit by a factor of 100, for example, 
barely affects the derived mass loss rate and changes the 
kinetic energy by an amount comparable to the uncer- 
tainty resulting from the power law fit in Eqn. [2] The 
general conclusion from Figure [2] is that, for values of 
a ~ 2-2.5, the GME mass loss rate for a saturated 
solar-type star is M ~ 5 x 10~^°Mq yr~^. The cor- 
responding CME kinetic energy requirement approaches 
Eke ~ O.IL©. In the context of current ideas concerning 
mass loss and efficiency of magnetic energy dissipation 
on active late-type stars these values are extremely high 
and their implications are discussed in ^^i^ and ^4.31 
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The implied CME mass and energy loss rates in stars 
with magnetic activity significantly below the saturation 
threshold depend much more heavily on the X-ray band- 
passes and the luminosity that enters the normalisation 
factor in Eqns. [5]and|ni We use a relation between coro- 
nal temperature and 0.1-10 k eV X-ray lurainosity for 
solar-like stars based on those of lGuedel et al.l (|1997| ) and 
iTelleschi et all (|2005l) . 

= 6 X lO^^T'^-S erg s-\ (10) 

where the constant and power law index have been tai- 
lored slightly so as to represent the isothermal plasma 
temperature, T, that reproduces the observed hardness 
^■atios L(o.2°ri.ofctv) derived by ITelleschi et all (l2005h . 
The hard band adopted by those authors is similar to the 
GOES band (1.54-12.4 keV) and the relation in Eqn.jlOj 
provides a reasonably accurate means for scaling the 
broad-band X-ray luminosity to this harder bandpass. 
Model hardness ratios and bandpass conversion factors as 
a function of temperature were derived using the APEC 
radiative loss model, as implemented in PIMM^. The 
resulting scaling factor depends approximately linearly 
on logLx, varying from ~ 0.3 at Lx = 3 x 10 erg s^^ 
to 10-3 at Lx = 10^7 erg s-^. 

Using the Lx to GOES bandpass scaling factor, we 
can obtain the CME mass and energy loss rates as a 
function of broad-band X-ray luminosity. These are il- 
lustrated in Figure 131 Here, the minimum and maximum 
flare energies were assumed to be Emin — ^O^^Emax 
and Ejnax = lO'^ix, though again the results depend 
only weakly on the exact limits of integration. The CME 
mass loss rate as a function of X-ray luminosity can be 
approximated by a simple power law at higher energies 
and a polynomial over a larger energy range: 

logMc = -54.6 4- 1.48 log ; Lx > 10^^ erg s"^ 

log Mc ^ -1339+131 logLx-4.371og^ Lx +0.049 log^ 

where the latter third order relation should be valid for 
Lx down to 10^® erg s"^. 

4. DISCUSSION 

The reader familiar with the literature touched upon 
in fjljon winds from late- type main-sequence stars might 
view the mass loss rates derived in fj3| with incredulity. 
The CME kinetic energy loss rate approaching 1/lOth of 
the stellar luminosity also seems implausibly high when 
compared with radiative losses through X-rays for satu- 
rated stars. For a range of spectral types this saturation 
level is consistently close to lO-^Lbo;, with a scatter of 
a factor of 2-3 (e.g. lWright et al.ll201lh . The derivations 
themselves in S|3| are, however, completely straightfor- 
ward and the origin of the numbers is easy to under- 
stand simply though inspection of the solar CME data. 
At face value, uncertainties in the analysis resulting from 
the power law fitting are also fairly small, resulting in 
mass and kinetic energy loss uncertainties of factors of a 
few. However, it should be kept in mind that there are 
relatively few CMEs at the higher energy and mass end 
of the observed distribution and systematic uncertainties 
resulting from this likely dominate. 

^ http:/ /cxc. harvard.edu/toolkit/pimms.jsp 




1 Q-^' L 

10=^ 10^" 10=5 10- 

X-ray Luminosity (erg/s) 




^q27 .|q28 iq29 ^q30 



X-ray Luminosity (erg/s) 

Fig. 3. — CME mass (top) and kinetic energy (bottom) loss rates 
vs. broad-band X-ray luminosity. The solid curve represents the 
power law index a = 2.25. The grey shaded areas represent the 
range in the loss rates corresponding to the power law index range 
1.5 < a < 3.0. 

4.1. Mass loss 

Also shown in Figure [T] is the vector corresponding to 
a constant ratio of mass loss to GOES X-ray energy loss 
converted to loss rates, M = 10-^°{Lx/10^°)Mq yr-\ 
and the mean CME ejected mass vs. flare X-ray fiuence 
weighted by the CME-flare association fraction. For flare 
energies < 3.5 x 10^^ erg, the latter lies remarkably par- 
allel to the former (probably by coincidence, but if not 
it raises an interesting issue), and offset by a factor of 
5 or so. One can see that this translates directly to the 
derived mass loss rate for a flare-dominated corona of a 
few 1O-^°M0 yr"^, as found in fJ3]and shown in FigurejH 
with little dependence on the power law index a. 

The value of a essentially controls the weighting be- 
tween flares of lower or higher energy. This is evident 
from Eqn. [51 for a < 2, the X-ray luminosity tends to 
be dominated by the larger flares, whereas for a > 2 
smaller flares contribute the largest fraction to the ob- 
served emission. Since the flare association-weighted 
mean CME mass follows very closely a constant E vs. 
mass, the total CME mass loss depends only weakly on 
a. Above the energy at which the CME-flare association 
is unity (3.5 x 10^^ erg), the slope of the £^-mass relation 
is more shallow and we expect the total derived mass loss 
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to be lower for a < 2, as is borne out in Figure [H 

It is of interest to see how the CME scahng might 
apply to solar levels of activity. From Figure |3l we 
find that, for a fairly active solar X-ray luminosity of 
Lx ^ 10^^ erg s~^, the mass loss rate is Mc 4 x 
1O~^^M0 yr~^ — a few percent of the total solar mass 
loss rate, in bro ad agreement with the average CME mass 
flux assessed bv lVourlidas et al] (|2010t) between 1999 and 
2003. 

The conclusion from the mean E vs. mass relation for 
solar flares and associated CMEs is that, if active stars 
are dominated by solar-like flares and we can extrapolate 
the relation in Figure [1] to higher flare energies, CME 
mass loss rates are very large — up to four orders of mag- 
nitude greater than the present-day solar wind and scal- 
ing with 1-8 A band X-ray luminosity roughly accord- 
ing to M - 5 X 10"1"(Lx/103°)Mq yr'^. lAarnio et all 
(|2012l ) find similar numbers for T Tauri stars (7 x 10"^-^- 
2 X 10~^Mq yr~^) using a similar method to that em- 
ployed here. 

Can such high mass loss rates possibly be correct? 
For the most active stars, the CME energy requirements 
clearly pose a problem; we return to this in t J4.3l In 
fJT] we cited evidence that upper limits to wind-driven 
mass loss in active solar-type stars based on attempts 
to detect free-free ra dio emission were of the order of 
a few lO-"Af0 yr"! IGaidos et all (pOOOl) . This anal- 
ysis assumed a spherically-symmetric wind, but a su- 
perposition of many CMEs should produce a qualita- 
tively similar, though probably more clumpy and tur- 
bulent, outflow (see also ii4.2[) . The most active star in 
the sample was tt^ UMa, with M < 5 x lO'^M© yr~^ 
This star has a broad-band X-ray luminosity of Lx ~ 
1029 erg s-i (|Drake et al.lll994t ITelleschi et al.ll2005l) — 
an order of magnitude or so below the saturation level. 
Based on our Mc vs. Lx relation in Figure[31 this implies 
M ~ 3 X IQ-'^'^Mq yr~i— well within the (jCaidos et al.l 
12000^) upper limit. This is also consistent with the max- 
imum allowed mass l oss rate of < l O'^Afn for a 10^ K 
"wind" estimated by iLim fc Whiti (11996) based on the 
requirement of radio transparency consistent with radio 
d etections of active s tars. 

IWood et al.l (pOOlh and Wood et al.l (HOGS') estimated 
mass loss rates for a handful of stars of different activ- 
ity level using astrospheric Lya absorption. While their 
analysis assumed a spherically-symmetric outflow with 
a wind speed of 400 km s~^, we would again expect 
a similar observational signature from a superposition 
of CMEs. The Wood et al. measurements scale with 
wind ram pressure, so an outflow comprised of generally 
faster CMEs would imply a proportionately lower mass 
loss rate. The mass-we i ghted mean CME speed in the 
lYashiro fc Gopalswamvl ()2009() sample is 1015 km s~^, 
implying a lower mass loss by only a factor of ~ 2 . For a 
small sample of G and K dwarfs lWood et all (|2005[ ) found 
M (X i^^-34±o.i8^ ^^Qj.^ jg ^j^g surface X-ray flux. This 
is similar to the CME mass loss power-law relation we 
found in Eqn. [TT] While they caution against extrap- 
olating the relation to higher activity levels, scaling to 
the X-ray surface flux for a solar-like star correspond- 
ing to = 3 X W° yields M - IQ'^^Mq yr'^, only 
a factor of a few lower than our CME-based estimate. 



At a flux level of tt^ UMa, the relation corresponds to 
M ^ lO^^^Af0 yr~^ — again similar to the CME scaling 
since the power law relations are also very similar. 

The scaling of the solar wind to higher magnetic activ- 
ity levels remains a very uncertain endeavour owing to 
the lack of a comp r ehens ive theory explaining the solar 
wind itself. ICohenI (poTTf ) argues that scaling mass loss 
according to X-ray luminosity is misleading because the 
latter is dictated by the closed magnetic field, while the 
former is dominated by open flux. He argues that mass 
loss rate s are unlikely to exceed IO^^^Mq yr~^. The 
model of iCranmer k, Saail ()2011|) is driven by the energy 
flux of magnetohydrodynamic turbulence from the con- 
vection zone and the filling factor of open field. Their 
mass loss rate for a saturated solar-like star is a few 
1O~^^M0 yr~^, which is two orders of magnitude lower 
than the direct CME scaling. This suggests that the 
indirect wind observations of Wood et al. could in fact 
be observations of quasi-continuous CME mass loss for 
the more active stars of the sample, rather than a direct 
analogy to the solar wind. 

Taken at face value, the solar CME data, combined 
with currently scant data on stellar winds and models of 
wind mass loss, suggest that at the highest activity lev- 
els mass loss could be dominated by CMEs, with a grad- 
ual transition to wind-dominated mass lo ss toward lowe r 
activities. In the context of the study of ICohenI (j2011l) . 
unlike a steady wind, CME mass loss is expected to scale 
with X-ray luminosity because flares and CMEs gener- 
ally originate from active regions that are dominated by 
closed field. 

In the context of the mass budget for individual CMEs, 
Equation [T] implies that for stellar flares with a total 1- 
8 A fiare X-ray fiuence of 10'^'* erg, the mean ejected mass 
is about 4 x 10^^ g. If the CME source plasma resides in 
the corona, this mass is uncomfortably large. It corre- 
sponds to the entire mass of a corona with a scale height 
of O.li?0, an emission measure of a few 10^^ cm~3 (like 
that of 47 Cas B from analysis of ITelleschi et al.l 120051) 
and a quiescent active region-like density of 10"'^'' cm~'^. A 
fiuence of 10'^^ erg, even limited to the 1-8 A GOES band, 
is still quite a modest flare compared with the largest 
flares seen on the most active stars and on T Tauri stars, 
whose broad-band X -ray fluences can reach 10^'' erg (e.g. 
ISchriiver et al.ll2012l ). This corresponds to ~ 3 x 10^^ erg 
in the 1-8 A band based on the scaling derived in SJ31 and 
would imply a mean ejected mass a factor of 30 higher 
still. These large mass requirements suggest that the so- 
lar CME X-ray fluence-mass relation must break toward 
the largest flares. 

4.2. Early Faint Sun Paradox 

Could the mass loss through CMEs on an early ac- 
tive Sun be rele v ant to the "early faint Sun paradox"? 
iSagan &: Mulled ()1972D pointed out that the lower so- 
lar luminosity predicted by stellar evolutionary the- 
ory earlier in the history of the solar system implies 
that for contemporary albedos and atmospheric compo- 
sition global mean temperatures would have been be- 
low the freezing point of seawater until about 2.3 Gyr 
ago, in contradiction with geological evidence for liq- 
uid oceans. Possible solutions to this paradox in- 
clude higher concentrations of greenhouse gases and 
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aerosols (e.g. iSagan fc Mulled 119721: iKastind fTOOl . a 
lower glo bal albedo, e ither through less cloud cover- 
age (e .g. iShavivl I2003D o r a smaller continental land 
mass osing et al.l l2010[ ). An alternative solution is 
an early Sun more massive by several percent t hat has 



since been whittled down by rnass loss (e.g. iGuzik et al 



20071) 



1987HSackmann k Boothrovdll2003l: iMinton fc Malhotra 



IWood et al.l (|2002f ) note that their inferred steady wind 
mass loss rates are insufficient when combined with re- 
lations for the secular decline of X-ray surface flux: the 
cumulative mass loss from an age of 1 Gyr or so is much 
less than 1% . As noted in i i4.1[ we find a relation be- 
tween mass loss from CMEs and X-ray flu x consistent 
with the observed steady wind relation of iWood et al.l 
(|2002l ) — M oc L]^^ for active stars, compared with their 
oc ^^-S-^ito.is ^ solar- like star such as Get with 

an age of about 0.75 Gyr and Lx ^ 10^^ erg s^^, the 
mass loss rate is similar to that of tt^ UMa considered in 
gm M ~ 3 X lO-^^M© yr-^ Even if such a rate lasted 
for over a Gyr, it would amount to less than 0.01 A/©, or 
an order of magnitude less than required to resolve the 
early faint Sun paradox unilaterally. 

4.3. Energy loss, dynamo saturation and CME 
implications 

A key question begging from the beginning of 52] is 
what fraction of the stellar bolometric luminosity can 
be scavenged by magnetic processes that give rise to 
flares and GMEs? Saturation of magnetic activity for 
the most active stars pegs broad-band X-ray radiative 
losses at ^ lO^^ibo;, which has generally been used as 
a saturation energy dissipation rate for convection zone 
dynamos. However, X-rays represent only one aspect of 
the energy budget of solar flares. The total radiative and 
non-thermal energy can be factors of 10-100 higher than 
the GOES 1-8 A fluence, and possi bly as much as as 
sociated GM E kinetic energies (e.g. Woods et al.' '200 
EmsUe et al.l [2005 : Ravmond 2008; Krctzschmar et al, 
20101 lKretzschmaijr2011 '). This translates to factors of 
~ 3 — 30 higher than broad-band X-ray fluence for flare- 
like temperatures. If active stellar coronae are dominated 
by flares and their behavior is similar on active stars, 
then energy requirements of flares alone could amount 
to 1% or more of the bolometric luminosity. 

We found in ij3]that scaling solar GME kinetic energy 
to a flare-dominated corona at saturated activity level 
would require a staggering fraction of the stellar luminos- 
ity, approaching 10%. If we declare 10% of Lboi too high 
a fraction of the total stellar energy budget to expend on 
GMEs, the implication is that the solar GME data cannot 
be extrapolated to significantly higher energies in the way 
we have done in [JH Bot h GME speed and mass increase 
with X-ray fluence in the lYashiro fc Gopalswamvl ()2009f ) 
sample. To avoid energy budget problems, either the 
GME kinetic energy vs. X-ray fluence must flatten out 
toward higher flare energies, implying that GME speed, 
mass, or both flatten out, or the GME-flare association 
rate must drop back significantly below 1. It is commonly 
argued that flares without assoc iated GMEs are confine d 
by overlying magnetic field (e.g. iSvestka &: GliveHll992l ). 
Active regions on active stars could confine more ener- 
getic GMEs associated with stronger flares because of 



stronger magnetic fields or different magnetic topology. 

A fit to the mean GME speed, analogous to those in 
Eqns. [U and m finds vdE) = 3.6 x lO'^E^ '^'^ km s"!, 
and for the same 10'^'' erg flare fluence corresponds to 
11 ,000 km s~^. The highest spe ed in the GME sample 
of lYashiro fc Gopalswarn^ ()2009D is about 3000 km s^^. 
The kinetic energy problem would be largely resolved 
were the ejection velocity to level out at this value toward 
higher flare energies. 

The kinetic energy problem for very large GMEs is 
possibly related to the observed cutoff in solar ener- 
getic particle (SEP) fluence for particle energies above 
10 Me V first inferred from cosmogenic radionuclides 
bv ILing cnfelte r fc Hudso n ( 1980, see also Ee cdy 199^ 
Hudsonll2007tlSchriiver et al.ii2012HUsoskin fc Kovahsovl 



20121) . The SEP fluence frequency spectrum breaks at 
approximately 10^*^ cm~^. Several different explanations 
have been suggested for this, such as an event energy 
dependence of SEP spectral distributions, particle prop- 
agation effects in the heliosphere, or SEP opening an- 
gles depending on the energy of the triggering event (e.g. 
iSchriiver et"al1 12012| ). [Hudsoni (|2007t ) notes that the 
SEP cutoff energy corresponds approximately to flares 
of XIO class (peak GOES 1-8 A flux of lO'^ Wm-^ 
or 3 X 10^^ erg s~^). He argues that there appears 
to be no corresponding cutoff in the frequency of such 
flares, and any steeping of the flare frequency spec- 
trum must happen at significantly higher energies. This 
suggests there might be another limiting factor govern- 
ing SEP production. Since SEPs are thought to be 
largely accelerated in GME-driven coronal and inter- 
planetary shocks (e .g., Kahl er et al.lll978l: lReameslll999l : 
iGopalswamv et al.l [2002: Oliver et al.ll2004[ ). a cutoff in 
energy could result from a corresponding limit to GME 
velocity or kinetic energy. How such a constraint might 
translate to significantly more active stars would then be 
important for the GME kinetic energy budget: whether 
it scales with maximum fiare, or active region, available 
energy, or is a more fundamental physical limitation with 
a cutoff at the same energy as seen on the Sun. 

We tentatively conclude that the relation between both 
GME mass and speed with flare fluence must flatten to- 
ward larger flare energies. This behavior would appear 
to differ to the scaling of magnetic flux a nd flare prop- 
erties from solar to active stellar cases. iPevtsov et al.l 
(|2003| ) find that, instead, the relationship between un- 
signed magnetic flux and X-ray spectral radiance for dif- 
ferent regions of the Sun scales over 12 orders of magni- 
tude to active stars. Flare temperatures, emission mea- 
sures and hard vs. soft components also appear to show a 
similar scaling fr om solar flares aU the way to giant flares 
on active stars (iFeldman et al.l [l995l : llsola et al.l [20071 : 
lAschwanden et al.ll2008[ )" Solar GME energies are lim- 
ited by the free energy available in solar active regions, 
which jGIopalswamy et al. (2010) note is < 10'^^ erg and 
usually at most lO^J ^-lO'^'^ erg (e.g. lEmsUe et afl 120121 : 
lAulanier et al.l [20131 ): in this context it would be of in- 
terest to assess the energy available in stellar active re- 
gions, which might be possible through Zeeman-Doppler 
imaging magnetograms. 

The large "hidden" GME and flare energy require- 
ments of a corona whose X-ray emission is dominated by 
flares implies that stars at the X-ray activity saturation 
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threshold are extracting niuch more than 10~^Li,oi of en- 
ergy, probably by an order of magnitude or more. This 
suggests that saturated stars are experiencing a magnetic 
energy dissipation limit, rather than a limit imposed by 
the ability of the star to sustain X-ray emitting loops due 
to centrifugal strip ping or poleward migration of mag- 
netic flux (StcpiciLeLaLl 120011: iJardine fc Unruhl 119991: 
iWright et~al.ii2011l) . That is, the maximum amount of 
the total energy budget able to be extracted by a mag- 
netic dynamo has been reached. If we allow CMEs to 
consume 1% of the stellar energy budget on X-ray satu- 
rated stars, the implied mass loss for saturated stars is 
Mc 5 X 10-"Mo yr-i. 



lations cannot be extrapolated to arbitrarily high flare 
energies: CME mass and/or kinetic energy vs. flare X- 
ray fluence must flatten off at flare energies > 10"^^ erg. 
A more reasonable CME energy budget of 1% of L^oi im- 
plies Mc ~ 5 X 1O~^^M0. Even for budgets an order of 
magnitude lower it seems likely that mass loss from active 
stars will be dominated by CMEs. The large "hidden" 
energy budget of flares and associated CMEs raises the 
question of what is the maximum amount of energy a 
solar-like star can extract from a magnetic dynamo? If 
saturated stars are dominated by flares, this energy is 
likely to be an order of magnitude larger than the ob- 
served broad-band X-ray saturation level of 10~^Li,oi. 



5. CONCLUSIONS 

Mean solar CME mass and kinetic energy are related 
to the associated flare X-ray fluence by power laws. If ac- 
tive stellar coronal X-ray emission is comprised of flares 
as observations suggest, and these flares adhere to the 
observed solar flare-CME scalings found here, very high 
CME mass loss rates exceeding 10~^°Mq are implied for 
the most active stars, consuming a tenth of the stellar 
bolometric luminosity. Since this energy requirement 
seems too high, we conclude that solar flare-CME re- 
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